Application of random phase approximation to vibrational excitations of double-A 

hypernuclei 
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Using the Hartree-Fock plus random-phase-approximation (HF+RPA), we study the impurity 
effect of A hyperon on the collective vibrational excitations of double-A hypernuclei. To this end, we 
employ a Skyrme-type AA'^ and AA interactions for the HF calculations, and the residual interactions 
for RPA derived with the same interactions. We find that inclusion of two A hyperons in ^^O 
shifts the energy of the collective states towards higher energies. In particular, the energy of the 
giant monopole resonance of aaO, as well as that of ^^Pb, becomes larger. This implies that the 
effective incompressibility modulus increases due to the impurity effect of A particle, if the /3-stability 
condition is not imposed. 

PACS numbers: 21.80.+a,21.60. Jz, 21.10.Re,23.20.-g 



I. INTRODUCTION 

Information on the interaction between a A hyperon 
and a nucleon deepens our understanding of baryon- 
baryon forces and the equation of state (EOS) of nuclear 
matter. In principle, the interaction between two parti- 
cles can be investigated with a measurement of their scat- 
tering. However, due to the short life-time of A hyperon, 
it has yet been difficult to perform a direct scattering 
experiment of nucleon and A hyperon. Therefore, the 
AiV interaction has been mainly investigated by 7 spec- 
troscopy of single- A hypernuclei [H-lll- Such measure- 
ments have revealed the A-impurity effect, that is, the 
change of several properties of atomic nuclei, such as ex- 
citation energies and transition probabilities of 7-ray, due 
to an addition of A particle. Apparently, high-resolution 
7-ray measurements are vital in investigating A hyper- 
nuclei. In addition to the existing experimental data, 
research projects currently planned at the J-PARC facil- 
ity using new Ge detector arrays (Hyperball- J) Q aim at 
obtaining new data on the low-lying energy level scheme 
of A hypernuclei in the sd shell region, that will lead to 
further understanding of the AA^ and AA interactions. 

Several theoretical calculations have been carried out 
to analyze the relation between low-lying energy levels 
of single- A hypernuclei and the AA^ interaction 043 ■ 
These calculations have not only contributed to iden- 
tification of energy level schemes of single-A hypernu- 
clei, but have also predicted the A-impurity effect on the 
structure of single-A hypernuclei, e.g., shrinkage of the 
radius of ^Li from ^LillO^], which was subsequently ob- 
served experimentally 

Besides single-A hypernuclei, double-A hypernuclei 
have also been studied both experimentally and theo- 
retically. Similar to the AA^ interaction, information on 
the AA interaction can be deduced from observation of 7- 
rays emitted from double-A hypernuclei. However, until 
now double-A hypernuclei have been produced only in an 
emulsion, and at present emitted 7-rays are difficult to 



detect experimentally with high precision. In addition, so 
far observed double-A hypernuclei in the emulsion have 
been limited to five cases(jYAHe and [Hi [HI); 

the experimental data have been scarce. Therefore the 
theoretical approaches make an important tool to assess 
the A impurity effect on the structure of double-A hyper- 
nuclei as well as appropriate selection of a target nucleus 
for future experiments. Theoretically, the double-A hy- 
pernuclei have been investigated within the frameworks 
of ab-initio few body model[17j, shell model|il8i] and clus- 
ter model[l^. However, these theoretical approaches de- 
mand a huge computational power, and they may be dif- 
ficult to apply to heavy hypernuclei. 

In order to study systematically the A-impurity ef- 
fect, from light to heavy nuclei, a Hartree-Fock (HF) 
plus random-phase-approximation (RPA) method pro- 
vides one of the most suitable tools. This approach 
has been applied to study vibrational excitations of nor- 
mal nuclei (without hyperons) throughout the nuclear 
chart, starting with a single energy functional applicable 
in the whole range of the nuclear chart. In particular, 
the RPA has been successfully applied to descriptions 
of giant resonances of atomic nuclei. See Refs.[2(]«J23'| 
for earlier applications of Tamm-Dancoff approximation 
to (K^,7r^) and (7r+,K+) reactions, and of an RPA- like 
model to single-particle spectra of single-A hypernuclei. 
So far, the mean field approach has been extended to A 
hypernuclei in order to study the ground state properties 
24I429I. the p otential energy surface in the deformation 
plane j30l432| and fission barrier heights ^Ij. Concern- 
ing excited states, the low-lying excited states of ^^Mg 
have recently been calculated with a 5-dimensional (5D) 
collective Bohr Hamiltonian on the basis of the Skyrme- 
Hartree-Fock method [34] (see also Ref.fis'l for a recent 
application of anti-symmetrized molecular dynamics to 
the A.^Mg hypernucleus) . Although the Bohr Hamilto- 
nian approach can handle a large amplitude collective 
motion, it is much easier to employ the RPA approach to 
describe collective vibrations with several multipolarities, 
including giant resonances as well. 
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In this paper, we extend the Skyrme-HF plus RPA 
(SHF+RPA) scheme to hypernuclei. Skyme-type AN and 
AA interactions, similarly to the Skyrme nucleon-nucleon 
interaction, are used in this work. The residual interac- 
tions for RPA are derived self-consistently from the sec- 
ond derivative of the energy functional with respect to 
densities. In this study, we shall focus on the double- 
A hypernuclei rather than single-A hypernuclei, partly 
because the description is much simpler due to the time- 
reversal symmetry. The A-impurity effect is expected to 
be stronger in double-A hypernuclei, and such calcula- 
tions will provide the upper limit of the impurity effect 
for single-A hypernuclei. 

The paper is organized as follows. In Sec. II, we de- 
scribe the formalism of the SHF+RPA for hypernuclei. 
In Sec. Ill, we apply the SHF+RPA method to ^^^O 
hypernuclei and present the results for the strength dis- 
tributions and the transition densities for the isoscalar gi- 
ant monopole resonance (IS 0+), the electric dipole (El), 
quadrupole (E2) and octupole (E3) transitions. We also 
calculate the isoscalar giant monopole resonance of j^^^Ph 
hypernucleus and discuss the nuclear incompressibility in 
the presence of A hyperon. We give a summary of the 
paper in Sec. IV. 



II. RPA FOR HYPERNUCLEI 

In order to describe the ground state and excited states 
of double-A hypernuclei, we adopt the Skyrme-type zero- 
range force for the AiV and AA interactions. The AA^ 
and 3-body ANN interactions of this type were first in- 
troduced by Rayet as[24], 

VAN{rA - Tat) = ^0 (1 + X^Pa)5{rA - Vn) 



S{rA - rN)k' 



t'^k'6{rA -rN)-k + iW^k'6{rA - ta,) • (cr x fc 



interaction, which originates mainly from the AA — 
coupling [35] . 

Together with the Skyrme NN interaction[36], the 
total energy Etot in the Hartree-Fock approximation is 
given by 



where 



Etot — Ejsf + Ea, 



En = j HNir)dr, 



(4) 



(5) 



is the energy for the core nucleus without A hyperons 
while 



-BA = y" [HNAir) + HAAir)]dr, 



(6) 



is due to the AA^ and AA interactions (see Appendix A) . 
The kinetic energy density for A particles is included in 
the energy density His[A{r). 

The SHF equations are obtained by taking variation 
of the total energy Etot with respect to the densities for 
neutrons, protons and A hyperons. These are given as 



-V- 



2TO*(r) 



V + UbN{r) + UbAir) Mr) = ebM^) 



(7) 

where the index b refers to proton, neutron or A, and 
£(, is the single-particle energy. The explicit forms for 
the mean-field potentials J7bjv(r) and C/{,a('"): and the 
effective mass TO^(r) are given in Appendix A. 

After we construct the ground state in the Hartree- 
Fock approximation, we describe excited states with RPA 
as a linear superposition of 1 particle-1 hole (Iplh) con- 
figurations. That is, the excitation operator Q^, for the 
fc-th RPA phonon is assumed to be. 



(1) 



i: ( 



vik) t 



(8) 



and 

VANN{rA,rNi,rN2) = 4^(''A - 'rNi)5{rA - t-atJ, (2) 

respectively. In a similar way, Lanskoy introduced the 
AA interaction as 129, 



waa(t'Ai - fAs) = Ao(5(rAi 

1/2 

+ -Ai [k' S{rAi - rA, 

+ A2fc'i5(rAi - taJ • fc 

'tAi + rA 



+ A3 



Pn 



- S{rAi - rAjfc' 



Sir A, - taJ. 



(3) 



The operator k' = — (Vi — V2)/2i acts on the left hand 
side while fc = ( Vi — V2)/2i acts on the right hand side. 
PN{r) is the density distribution for the nucleons. The 
last term in Eq.(I3]) corresponds to the three-body AAA^ 



where X^'^ and Y^^"^ are the forward and backward am- 
plitudes, respectively. a|, and a\ are the creation opera- 
tors for a particle state p and for a hole state h, respec- 
tively. The excitation energy Ek is obtained by diagonal- 
izing the 2i^-dimensional RPA equation. 



A B 

-B* -A* 



(9) 



where v is the number of Iplh configurations. Here, A 
and B are RPA matrices given by. 



Aph^p' h' 
^ph^p' h' 



y^p 
^pp' hh' 



£h)Spp'6hh' + Vph'hp' 



(10) 



where v is the residual interaction derived from the en- 
ergy functional, -Etot- The formalism is almost the same 
as the standard RPA found e.g., in Refs.^371. Issj . but 
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the particle-hole configurations run over not only pro- 
tons and neutron but also A hyperons. The interaction 
matrix elements Vph'hp' and Vppihh' include the A7V and 
AA interactions as well as the NN interaction (see Ap- 
pendix B). 

The external fields for electric multipole excitations 
with multipolarities L =i and 1 are defined as 



Fel = e^^rf Ylm (n), 



(11) 



while that for the "isoscalar" monopole transition is 

(12) 



For the electric dipole response, we take into account the 
center of mass motion of the whole hypernucleus and use 
the operator 



Fei 



eY,{nYui{h) - RYim{R)), 
NniN + Natua 




ieA / 



where 



H= — 

M 



(13) 



(14) 



is the center of mass of the hypernucleus, and M = 
mAsiZ + A^) -I- uiaNa is the total mass, tun = {rrip + 
m„)/2 = 938.92 MeV/c^ and tua = 1115.68 MeV/c^ be- 
ing the mass of nucleon and A hyperon, respectively. N, 
Z and Na are the number of neutron, proton and A hy- 
peron, respectively. 



III. RESULTS 

A. single-particle level of aaO 

We now numerically solve the RPA equation and dis- 
cuss the collective excitations of double-A hypernuclei. 
Before we show the results for multipole vibrations, we 
first discuss the single particle levels of double-A hyper- 
nucleus j^^O and ^^O, which will help to understand the 
A-impurity effect on the giant resonances. To this end, 
we assume spherical symmetry, and solve the SHF equa- 
tion in the coordinate space with a grid size of dr=0.1 
fm. We use the SkM* parameter set for NN interac- 
tion |39!|, while the No. 5 parameter set in Ref.[25'] for the 
AA^ interaction, whose parameters were determined by 
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FIG. 1: (Color online) Neutron and proton single-particle lev- 
els of (the solid lines) and aaO (the dashed lines) obtained 
with the Skyrme-Hartree-Fock method. 



fitting the Hartree-Fock calculations to the experimental 
binding energies of single- A hypernuclei [25.]. For the AA 
interaction^ we use the SAAl parameter set evaluated by 
Lanskoy[26l|. This parameter set was obtained by fitting 
to the AA bond energy [26^], ABaa = Baa — 2-Ba, where 
Ba is the one- A hyperon separation energy from a ^'''^Z 
hypernucleus and -Baa is the two-A hyperon separation 
energy of ^aa"^- ^^^^ show in the next subsection, 

the dependence of giant resonances on a choice of param- 
eter sets for the AA^ and the AA interactions is weak, and 
any significant change in the results is not obtained even 
if we use different parameter sets for the interactions. 

Figure [T] shows the neutron and proton single-particle 
levels of the ^^O (the solid lines) and the aa^ (the dashed 
lines). The single-particle energies of the aaO hypernu- 
cleus are smaller than those of ^^O, since the depths of 
the central part of the mean-field potentials are deepened 
both for protons and neutrons due to the attractive AA^ 
interaction, as shown in Fig. 2 of Ref . (27j . The lowest 
level (lsi/2) is the most sensitive to the addition of A 
particles, for which the single particle levels are shifted 
by Ae„ = -1.4 MeV for neutron and Ae^ = -1.3 MeV 
for proton. This A-impurity effect becomes weaker as 
the energy of a single-particle state increases. If the con- 
tinuum spectra are discretized within a large box, the 
difference of single-particle energies between the and 
the AA^ nuclei is much smaller as compared to the bound 
levels. Consequently, in the independent-particle approx- 
imation, the excitation energies increase relatively in the 
double-A hypernucleus as compared to those of the nor- 
mal nucleus. In the next subsections, we will see that this 
is the case even in the presence of the residual particle- 
hole interaction. 
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TABLE I: The excitation energies and the electromagnetic 
transition probabihties, B{E2) and B{E3), for the first 2"*" 
and 3~ states of '^^O and aaO nuclei obtained with the 
Skyrme-HF+RPA method. 





2^ state 


3j state 


nucleus 


E (MeV) B{E2) (e^fm'') 


E (MeV) B{E3) (e^fm") 


ibQ 


13.1 0.726 


6.06 91.1 




13.8 0.529 


6.32 67.7 



B. Low- lying excitations 

We next solve the RPA equation in order to discuss col- 
lective excitations of the hypernucleus. To this end, 
we discretize the single-particle continuum states with 
the box boundary condition with the box size of 16.0 
fm. We take into account the continuum states up to 
e=30 MeV, and consider the Iplh configurations whose 
unperturbed energy, ep — eh-, is smaller than 60 MeV. For 
the residual interactions, we neglect the Coulomb and 
the spin-orbit terms for simplicity, although we include 
all the other terms self-consistently. Therefore, our RPA 
calculations are not fully self-consistent, and the spurious 
translational motion appears at a finite excitation energy. 
In order to recover effectively the self-consistency, we in- 
troduce a scaling factor / to the residual interaction Vres 
so as to produce the spurious translational mode at zero 
energy. 

Table |T] shows the results of such RPA calculations for 
the lowest quadrupole and the octupole states of ^^O and 
^^O. For both the 2^ and 3^ states, the impurity effects 
of A particles slightly reduces the collectivity, that is, 
the excitation energies are increased while the electro- 
magnetic transition probabilities are decreased by 26-28 
%. The increase of the excitation energies is consistent 
with the increase of unperturbed particle-hole energies 
discussed in the previous subsection. 



C. Giant resonances 

The RPA method has been successfully applied to gi- 
ant resonances of the normal nuclei. We therefore apply 
it in this subsection to the giant resonances of the double- 
A hypernucleus j^^O, although they may not be easy to 
access experimentally at present. The top, the middle, 
and the bottom panels of Fig. [5] show the strength dis- 
tributions for the electric dipole (El), quadrupole (E2) 
and octupole (E3) excitations, respectively, weighted by 
the Lorentzian function with a width of 1.0 MeV. The 
solid and the dashed lines denote the results for the j^^O 
and ^^O nuclei, respectively. In order to assess the role 
of A hyperon, we also show the results for the j^^O hy- 
pernuclei in which the AN and AA interactions are taken 
into account only in the ground state, that is, the results 
obtained by switching off the residual AA^ and AA in- 
teractions (the dotted lines). The figure indicates that 
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FIG. 2: (Color online) The strength distributions for the elec- 
tric dipole (El, the top panel), the electric quadrupole (E2, 
the middle panel) and the electric octupole (E3, the bottom 
panel) excitations of the ^^O nucleus (the dashed lines) and 
of the double-A hypernucleus aaO (the solid and the dotted 
lines). The solid lines are obtained by including the residual 
A'^A'^, AA'^ and AA interactions, while the dotted lines are ob- 
tained by including only the residual A'A'^ interactions. The 
strength distributions are weighted by the Lorentzian func- 
tion with a width of 1.0 MeV. For the peaks indicated by the 
arrows, the transition densities are shown in FigO 



the addition of the A hyperons shifts the peaks of the 
strength functions towards high energy for all the modes 
of excitations. This is similar to the results for the low- 
lying modes of excitations discussed in the previous sub- 
section, and is again caused mainly by the change of the 
single-particle energies. On the other hand, the difference 
between the solid and the dotted lines is relatively small, 
except for the low-lying dipole state at i? = 12.8 MeV. 
We have confirmed that the strength functions remain 
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TABLE II: The centroid energy E for the El, E2 and E3 
modes of excitations for ^^O and aaO nuclei. Those are given 
in units of MeV, and the results of both the unperturbed (HP) 
and the perturbed (RPA) calculations are shown. 5E denotes 
the difference of the centroid energies between aaO and ^^0. 
The values in the parentheses for the El mode are the results 
obtained by excluding the contribution of the low-lying peak 
aX E = 12.8 MeV. 







El E2 E3 


(HE) 


ibQ 

AA*-* 


13.76 25.57 26.53 
14.34 26.63 27.74 


5E 


+0.58 +1.06 +1.21 


(RPA) 


ibQ 

18 O 


19.92 19.55 22.32 
19.68 (20.95) 20.09 24.05 


5E 


-0.24 (+1.03) +0.54 +1.73 



almost the same, including the low-lying dipole state, 
even if other parameter sets for the AiV and AA interac- 
tions are employed. This suggests that the main effect 
of A particles on the collective vibrational excitations is 
indeed attributed to the change in the single-particle en- 
ergies, rather than the residual AA^ and AA interactions, 
although the low-lying dipole state may require a sepa- 
rate analysis (see Fig. 4 below). 

In order to see the A-impurity effect quantitatively, 
we show in Table |lT] the centroid energy defined &s E — 
mi /mo, where is fc-th energy- weighted sum-rule. 



ruk 



(15) 



for the unperturbed (HF) and the perturbed (RPA) 
strength functions. We also list the difference of the cen- 
troid energy, 5E, between the ^fO and the ^^O nuclei. 
The values in the parentheses for the RPA El response 
are the results obtained by excluding the contribution of 
the low-lying dipole peak at E — 12.8 MeV. As is ex- 
pected, the centroid energies E for the HF calculations 
shift to higher energies when the A hyperons are added to 
^^O. For the E2 and E3 responses, this remain the same 
even if the residual interactions are taken into account in 
RPA. For the El response, the energy shift is negative, 
but it turs to positive if the low-lying peak is excluded. 
As we will show below, this low-lying peak corresponds 
to the dipole motion of a A hyperon against the core 
nucleus. We thus conclude that the A hyperons gener- 
ally increases the centroid energy for collective motions 
of the core nucleus, not only for the quadrupole and the 
octupole states but also for the dipole states. 

Figure [3] shows the transition densities for the giant 
dipole and quadrupole resonances and for the high-lying 
octupole state, which are indicated by the arrows in 
FiglH The top, the middle and the bottom panels denote 
the transition densities for the neutrons, the protons and 
the A hyperons, respectively. Those transition densities 



are computed as 

ph 

X ipp{r)Lph{r){jplp\\YL\\jhlh), 



(16) 



where j^, li are the total and the orbital angular mo- 
menta for a single-particle state «, respectively, while 
ipi (r) is the radial part of the wave function, normalized 
as ip*{r)tpi{r) r'^dr = 1. The peaks of the transition 
densities for the neutrons and protons slightly move to 
a small distance and are shifted towards inside for all 
the multipolarities when two A hyperons are added. The 
amplitude of the transition density for the A hyperons is 
about 1/10 — 1/100 smaller than that for the protons and 
neutrons, so that the A hyperons do not contribute much 
to these giant resonances. For the E2 and E3 states, the 
neutrons and the protons oscillate in phase as is expected 
for isoscalar motions, while they oscillate out of phase for 
the El state [i.e., the isovector motion). In addition, the 
A hyperons oscillate in phase with the protons and the 
neutrons for the E2 and the E3 modes, while they oscil- 
late in phase with the protons for the El mode. When 
the Coulomb force is turned off completely, that is, when 
the single-particle levels for the protons are identical to 
those for the neutrons, the amplitude of the El transition 
density for the A hyperons vanish. 

The low-lying dipole state at E — 12.8 MeV deserves a 
special attention. This peak appears only when the A hy- 
perons are added to the ^^O nucleus, and a similar peak 
does not seen in other modes of excitations. Figure |4] 
shows the transition density for this state. In contrast to 
the giant resonance shown in Figl21 the amplitude of the 
transition density for the A hyperons is higher than that 
for the protons and the neutrons. The strongest RPA 
amplitude, ^ = — Y^, contributing to this peak is the 
[lp(ls)~^] configuration of the A particles (f — 0.873). 
The total RPA amplitudes for the neutrons and the pro- 
tons are small (^ = 0.050 for the neutrons and ^ — 0.071 
for the protons), and these values become entirely zero 
when the AN interaction is switched off. The neutrons 
and the protons oscillate in phase, and the A particles 
oscillate out of phase with the nucleons. We can thus 
interpret this mode as a dipole oscillation of the A par- 
ticles against the core nucleus ^^O, similar to the soft 
dipole motion of a valence neutron in halo nuclei (40| . 



D. Giant monopole resonance and 
incompressibility 

Giant monopole resonances, the so called "breathing 
mode", are intimately related to the incompressibility of 
nuclear matter |4ll447| , which plays an important role in 
neutron stars. It has been shown that the EOS of infi- 
nite nuclear matter is softened when hyperons(A, S, E) 
emerge at high densities, and as a consequence the max- 
imum mass of neutron stars becomes smaller [48145 Oll . It 
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FIG. 3: (Color online) The transition densities for the giant dipole (the left panel) and the giant quadrupole (the middle panel) 
resonances and for the high-lying octupole state (the right panel) in aaO (the solid lines) and '^^O (the dashed lines). The 
corresponding states are denoted by the arrows in Fig. [51 The top, the middle, and the bottom panels denote the transition 
densities for the neutrons, the protons and the A hyperons, respectively. 
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FIG. 4: (Color online) The transition density for the low-lying 
El state at E = 12.8 MeV in the /aO hypernucleus. The 
dashed, the dotted, and the solid lines show the contributions 
of the neutrons, the protons, and the A hyperons, respectively. 



is thus of interest to investigate the A-impurity effect of 
giant monopole resonances in finite nuclei. 

Figure \5\ shows the strength function for the isoscalar 
monopole responses of ^^O and ^aO nuclei. The meaning 
of each line is the same as in Fig. [21 For a comparison, 
the figure also shows the strength function for ^"^Pb and 
As in the other multipolarities discussed in the 



previous subsection, the strength distributions are shifted 
towards high energies when A hyperons are added, and 
also the difference between the solid and the dotted lines 
is small, indicating that the residual AA^ and AA inter- 
actions play a minor role. 

Figure [Bl shows the transition densities for the giant 
monopole resonances corresponding to the states indi- 
cated by the arrows in the Fig [51 (that is, those states 
at £; = 22.2 MeV and 21.1 MeV for ^^^O and ^^O, re- 
spectively, and at 14.6 MeV and 14.2 MeV for f^Ph and 
^•^^Pb, respectively). The meaning of each line is the 
same as in Fig. [31 For the oxygen nuclei, when A hyper- 
ons are added, the amplitude of the transition density 
for the neutrons decreases by about 20% while that for 
the protons remains almost the same. The amplitude of 
the A transition density is about 10 times smaller than 
that of the nucleons. It is interesting to notice that the A 
hyperons oscillate out of phase with the nucleons. These 
features are qualitatively the same for the lead nuclei as 
well, although the changes in the transition densities are 
much smaller compared to the oxygen isotopes. 

According to Blaizot, the effective incompressibility 
modulus Ka for finite nuclei without A hyperons is de- 
fined as liTI. 
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FIG. 5: (Color online) The strength distributions for the 
isoscalar monopole mode for the and aaO nuclei (the 
top panel) and for the ^"^Pb and ijvPb nuclei (the bottom 
panel). The meaning of each line is the same as in Fig. [51 



TABLE III: Properties of the isoscalar monopole responses 
obtained with the Skyrme HF+RPA method. iJo+ is the cen- 
troid energy, and £ is defined as £ = -\/mi/m_i. (r^)n 



+p 



and a/ (r^) are the root mean square radii for the core nu- 
clei and that for the total densities, respectively. Ka is the 
effective nuclear incompressibility defined by Eqs. (I17II and 
HI. 







^(r2)(fm) 


A'A(MeV) 


ibQ 


22.4 


21.7 


2.68 


2.68 


81.6 




24.3 


23.5 


2.64 


2.58 


90.0 


2U8pb 


14.1 


14.0 


5.56 


5.56 


146 




14.5 


UA 


5.55 


5.53 


153 



^) is the root mean square radius, and £^ 



where • 

mi/m_i (see Eq. (|15|) for the definition of the k-th 
energy-weighted sum-rule, m^). When the A hyperons 
are present, this formula is modified as. 



Ka 



N + Z{r^)n+p mNNA{r^)A 



A 



(18) 

where we have used Eqs. (3.45) and (3.47) in Ref.|4l| and 
the energy-weighted sum-rule for the isoscalar monopole 
transition, 
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FIG. 6: (Color online) The transition densities for the 
isoscalar giant monopole resonances of aaO (the left panels) 
and AyvPb (the right panels), with comparisons to the transi- 



tion densities for and 
the same as in Fig. |3l 



208 



Pb. The meaning of each line is 



mi(L = 0) = {N + Z)(r2)„+p + —NA{r')A. (19) 

mpf niA 



In Eq. 



is the root mean square radius 



of the core nucleus. Notice that Eq. (|18l) is reduced to 
Eq.inD when_7VA = 0. In Table Hill we hst the cen- 
troid energy Eq+ , £ for the isoscalar monopole modes, 
the root-mean-square radii, (r^) and -y/ {r'^)n^ 



and 



the effective incompressibility, Ka, calculated according 
to Eqs. pT|) and ^TE\i . When A hyperons are added, the 
centroid energies increase by 1.9 MeV for ^^O and 0.4 
MeV for ^"^Pb, and the rms radii for the core nucleus, 



/ii+pj 



decrease by 0.04 fm for ^^O and 0.01 fm for 



^"^Pb. As we have shown, the increase of the centroid 
energies is mainly due to the change of single-particle lev- 
els, and the residual AA^ and AA interactions give only a 
minor effect. The decrease of the rms radii is attributed 
to the attractive AN interaction, that is, the shrinkage 
effect of A hyperons. The effective incompressibility, Ka, 
increases for both the nuclei studied when A hyperons are 
added. 

The increase of the effective incompressibility should 
reflect the properties of infinite nuclear matter. In or- 
der to assess this. Fig. 7 shows the binding energy per 
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30 




P (fm"') 



FIG. 7: (Color online) The binding energy per particle in 
infinite nuclear matter for the fraction of A particle of a;A=0 
(the dashed line) and a::A=0.11 (the solid line). The neutron 
and the proton densities are set to be equal, that is, pn ~ 
Pp = (1 - xa)p/2. 



particle in infinite nuclear matter, E/A, given as 



E 
A 



H 
P 



P 



2mA 



1 



PNPA 



Ao 



1 



, ^3 2 a 

+^PaPn 



3X2)pata 



(20) 



where Hn is the nucleon part of the energy density eval- 
uated in infinite matter, p is the total density, while 
Pa = x\p and pN = (1 — x\)p are the densities of A 
particles and nucleons, respectively. The kinetic energy 
densities rjv and ta are evaluated as tn = •^PNkpj^ /5 and 
ta = Spa^fa/^j respectively, where the Fermi momenta 
are given by fc_F AT = (37r^piv/2)^/^ and kpA — (Stt^pa)"^^^- 
The dashed and the solid lines in Fig. 7 show the results 
with the A fraction of a;A=0 and a;A=2/18, respectively. 
Here, we have assumed that the neutron and the pro- 
ton densities are the same, Pn = Pp = Pn/'2- One can 
see that the addition of A particles shifts the equilibrium 
density po towards a high density, that is, po = 0.161 
fm-3 for XA=0 and po = 0.185 fm^^ for xa=2/18. The 
incompressibility for infinite nuclear matter. 



, {<P{E/A)\ 



dp' )o= 



(21) 



P=PO 



is = 217 MeV for xa=0 and K^o = 239 MeV for 
x\ — 2/18, agreeing with the increase of effective incom- 
pressibility of the finite hypernuclei. 

Notice that this observation does not contradict to the 
fact that hyperons soften the equation of state for infi- 
nite nuclear matter relevant to neutron stars. That is. 



in neutron star calculations, the emergence of hyperons 
takes place at high densities and nucleons are the only 
constituents at the normal density, when the beta sta- 
bility condition is imposed. In contrast. Fig. 7 shows 
the effect of hyperons on the incompressibility defined at 
the equilibrium density. Even though the beta stability 
condition does not hold there if the A fraction is finite, 
this EOS is more relevant to giant monopole resonances 
of finite hypernuclei. 



IV. SUMMARY 

We have extended the Skyrme-HF plus RPA schemes 
to calculations for vibrational excitations of double-A hy- 
pernuclei. We have applied it to the electric dipole (El), 
the quadrupole (E2), and the octupole (E3) modes of 
excitations in the aa^^ hypernucleus. We have shown 
that the strength distributions shift towards high ener- 
gies for all the modes when the A hyperons are added to 
^^O. This is the case both for the low-lying quadrupole 
and octupole states and for giant resonances. At the 
same time, the electromagnetic transition probabilities 
also decrease. We have argued that these features are 
mainly caused by the change in the single-particle ener- 
gies, whereas the residual AiV and AA interactions play a 
minor role. The calculated transition densities show that 
the peak of the transition densities are shifted towards 
inside and the height of the peaks slightly changes due 
to the impurity effect of A hyperons. 

For the El strength, we have found a new peak at low- 
lying energy, that is absent in the El response of the core 
nucleus. From the analysis of the transition density, we 
have shown that this state corresponds to an oscillation 
of the A particles against the core nucleus. 

We have also discussed the A-impurity effect on the 
isoscalar monopole vibration of aa^ and aa Pt), and the 
incompressibility of infinite nuclear matter in the pres- 
ence of A hyperons. When the A hyperons are added, the 
strength distributions are shifted to higher energies and 
thus the centroid energies increase, similarly to the other 
multipole transitions. We have shown that the transi- 
tion density for the A particles behave rather differently 
from the transition densities for the neutrons and the pro- 
tons. The increase of the centroid energy for the giant 
monopole resonance implies that A particles increase the 
nuclear incompressibility, when hyperons were emerged 
at the equilibrium density. 

In this paper, we have studied several collective vibra- 



tional motions, taking the aa^ hypernucleus as exam- 
ples. It would be an interesting future work to study 
systematically the A-impurity effect on the collective ex- 
citations of other double-A hypernuclei. In particular, 
the low-lying dipole mode, originated from a dipole os- 
cillation of the A particles against the core nucleus, would 
be interesting to study. For this purpose, we would have 
to extend our formalism by including the pairing correla- 
tions with the quasi-particle RPA (QRPA). Another in- 



teresting extension is to study the collective excitations of 
single-A hypernuclei, although the broken time-reversal 
symmetry will have to be taken into account correctly 
there. 
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Appendix A: Energy densities and mean-field 
potentials for hypernuclei 

In this Appendix A, we summarize the explicit formu- 
lae for the energy densities -/Jata and i?AA in Eq.lIH]) and 
the mean- field potentials in Eq.©. 

The energy density H^k, due to the Skyrme-type A7V 
interaction given by Eq.((T]), reads (2^ . 



UAN{r) = 4 (l + PN + \ (t^ + 4) 



(A.4) 



and 



UAA{r) = ^AoPA + ^(Ai + 3A2)ta 



(A.5) 



Note that the index q refers only to the proton and the 
neutron. The effective mass for the nucleons and the A 
hyperons in Eq. ([7]) are given by 

1 



2m* 2mAr 4 



+ 7(^1 +*2)paW, (A.6) 



and 



2m\ 2mA 4 
respectively. 



+ -M + t2)PN{r) + ^(Ai 3A3)pA(r), 



(A.7) 



HNA{r) = TA + to 1 + o^O PNPA 

2mA \ 2 J 

+ 4){rAPN + TAT pa) + ^(3t^ - 4){VpN ■ Vpa) 

+ l^^C^PN ■ J A + VpA • Jn) + \taPK{pl{ + 2p„Pp), 

(A.l) 

while the A A part, Haa, originated from the A A interac- 
tion given by Eq.® reads |26|. 



HAA{r) = ^AoPa + -^(^1 + 3A2)pATA 
+ ^(Aa - Ai)paVVa + \xzpIp%- 



(A.2) 



Here, pf, = Pbir), n = Tf,(r), Ji, — Jb{r) are the num- 
ber, the kinetic energy and the spin densities, respec- 
tively (6 = P^n, or A). The indices p, n and A are 
the nucleon, the proton, the neutron and the A hyperon, 
respectively. 

After taking variation of the energy in Eq.(|3]) with re- 
spect to the densities, we obtain the Skyrme-Hartree- 
Fock equation given by Eq. ([7]) . The mean- field potentials 
in Eq.([7]) are given by 

U,A{r) =t^[l + PA + \ {t^ + 4) TA 

- I (3t^ - 4) VVa - ^< V . J A (A.3) 

+ ^W^o'^VpA • {-iV X cr) + ^4pa{2pn - Pq), 



Appendix B: AA'^ and AA residual interactions 

The matrix elements for a particle-hole residual inter- 
action Wrcs are given as(37l [38| 



vph'hp' - {pih)-^LK\v,,,\p'{h')-^LK), 
vpp'hh' = {pih)-^LK,p'{h')-^LK\vros\), 



(B.l) 



where L is the multipolarity for the particle-hole ex- 
citations and K is its z-component. For hypernuclei, 
the residual interaction can be separated into two parts. 



Vrc 



,,6if)2 



(N) 



,,fcl&2 



(A), where the indices bi and 62 



denote p, n or A. The interaction v^^^^{N) is due to 
the NN residual interaction, whose explicit form can be 
found in e.g. Refs.fsi'-'Hs'l. v^.^^^{A) is the additional term 
due to the AA^ and the A A residual interactions. These 
are given in the form of. 



-I- C6jb2(Vi - 



Vf + V^2-(Vi-V'i)(V2- V^) 



V'i)(V2+V'2) 



(B.2) 



where a^^bj, bi,^i,2 and Cb^b2 are given by 



PA (3 



<J ■ <J T ■ T 



— Cqq' = 0, 



(B.3) 
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+ 4 {pN-^)+^ap%-'pA 

1 



(B.4) 



for (g, g'' = p or n), and 

AAA = ^-^0(1 - (Tl ■ 0-2) + ^A3P^(1 - <Ti • 0-2) 
&AA = (Ai(1 - <7i • 0-2) + A2(3 + CTi • 0-2)) (B.5) 

- CTi ■ (T2) - 3A2(3 + <7i ■ (7-2) 



CAA 



for the A A terms. 
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